Observation of a One-Dimensional Spin-Orbit Gap in a Quantum Wire 
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Understanding the flow of spins in magnetic layered structures has enabled an increase in data 
storage density in hard drives over the past decade of more than two orders of magnitude 0]. 
Following this remarkable success, the field of 'spintronics' or spin-based electronics [l|, 2, H| is 
moving beyond effects based on local spin polarisation and is turning its attention to spin-orbit 
' interaction (SOI) effects, which hold promise for the production, detection and manipulation of 

\ spin currents, allowing coherent transmission of information within a device P, Q|. While SOI- 

induced spin transport effects have been observed in two- and three-dimensional samples, these 
have been subtle and elusive, often detected only indirectly in electrical transport or else with more 
sophisticated techniques 0, [H, [f| 0, HI, 01 ■ Here we present the first observation of a predicted 'spin- 
orbit gap' in a one-dimensional sample, where counter-propagating spins, constituting a spin current, 
■ are accompanied by a clear signal in the easily-measured linear conductance of the system [loL Hl|. 

<D 

We first introduce the class of phenomena we dub 'the one-dimensional spin-orbit gap' using a simple example 
adapted from Ref. [Io| , then describe our experiment in detail, and finally present a more elaborate model which 
. captures most of the features seen in our data. 

The spin-orbit interaction is a relativistic effect where a charged particle moving in an electric field experiences 
an effective magnetic field which couples to its spin In semiconductor heterostructures, the electric field can 

arise as a result of either the lack of an inversion centre in the crystal (bulk inversion asymmetry, BIA) or a lack 
of symmetry in an external confining potential (structural inversion asymmetry, SIA) due to crystal interfaces or 
q ' additional structures such as metallic gates [l3j] • The strength of the resulting effective magnetic field is proportional 
l — to both the particle's momentum and the original electric field. 

i Consider a spin-degenerate one-dimensional subband with a Hamiltonian Hq = , where h is Planck's constant, 
k the particle's momentum and m its mass (Figure la). The leading order SO contribution to the Hamiltonian is 

Hso = ■ { k x X7V), where a is the particle's spin, V the electrostatic potential and (5 a material-dependent 
parameter [14j • This term breaks the spin-degeneracy of the system and results in two spinful subbands separated by 
a lateral (wave vector) shift as shown in Figure lb. 

Despite this rather striking change in the band structure, measurements of conductance through the system cannot 
distinguish the situation shown in Figure lb from the case where the spins are degenerate. In both cases the edges 
of the two spin subbands occur at the same energy, so in both cases the conductance rises by Go = 2e 2 /h when the 
Fermi level of the system is tuned through this energy (e.g. by applying a voltage to a nearby gate) [HI, [3, 
(Figure Id.) To detect the SO shift in a transport measurement, a different approach is needed. 

Note that the spins as shown in Figure lb are polarized in the direction of B so, which is perpendicular to both 
the momentum (i.e. the ID wire) and the external electric field [HI]. A magnetic field applied exactly along B so 
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■ shifts the spinful bands up and down by the Zeeman energy respectively, splitting each step of size Go into two steps 
of size Go/2 as shown in Figure le. Here spin is polarized alternately up and down for the respective spinful bands 
and thus the charge current can be completely spin-polarized when only one spinful band is occupied; however, charge 
transport measurements still reveal no difference between this system and one where there is no SOI. 

Now consider a magnetic field applied perpendicular to B so (say along the wire): the two spinful bands are mixed 
by this term so that the zero field crossing-point becomes an anti-crossing. We call this feature the 'spin-orbit gap' 
(Figure lc). When the Fermi level lies within such a gap, two Fermi points (e.g. 3 and 4 in Figure If) contribute to 
the current as opposed to the four at the same energy in Figure Id. The conductance is thus reduced by Go/2 as 
shown in Figure 2d. In addition, the current from this subband is completely spin polarized in a way that is expected 
to be robust to moderate disorder: holes must scatter between the points 1 and 2 or between 3 and 4 (Figure If) 
and perform a spin flip in order to backscatter. In such a system of counter-propagating spins (very nearly what is 
termed a 'helical liquid' [19]), the direction of the spin current is expected to be independent of the sign of the voltage 
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FIG. 1: The spin-orbit gap in a simple model and the associated conductance features, a, The dispersion relation for electrons 
in a spin-degenerate one-dimensional subband. b, The spin-orbit interaction (SOI) lifts the spin degeneracy, displacing spinful 
subbands laterally with respect to each other, c, An applied magnetic field can mix the two spin bands, creating the anti- 
crossing we dub the 'spin-orbit gap'. When the Fermi energy is tuned to be within this gap, particles of opposite spin travel 
in opposite directions, producing a spin current. At the same time, a clear drop is expected in the conductance, shown in f. 
d, With SOI but without an applied magnetic field, the conductance of the system increases by a step of 2e 2 jh each time the 
gate voltage aligns the Fermi energy to the bottom of a subband. e, A magnetic field applied parallel to B so splits each step 
into two half-steps of e 2 /h. f, With a magnetic field applied perpendicular to Bso, spin-orbit gaps appear as in c and the 
conductance drops when the Fermi energy lies within such a gap. 

applied across the system and is determined only by the sign of the SOI. Intriguingly, a pure spin current, without 
charge current, is expected to exist at zero bias voltage. 

7- 

Finally, a magnetic field applied neither exactly perpendicular nor parallel to B so results in a mixture of the two 
effects described above [Io| . 

We turn now to our experiment. Our samples are GaAs/AlGaAs hole quantum wires produced by the cleaved 
edge overgrowth (CEO) method [2^, [2l[. Starting with an extremely high-mobility 2D hole gas (2DHG) realised in a 
carbon-doped AlGaAs/GaAs/AlGaAs quantum well [22|], the sample is cleaved and additional AlGaAs is grown using 
Molecular Beam Epitaxy (MBE) over the freshly exposed surface. Further carbon-doping leads to hole accumulation 
at the GaAs/AlGaAs interface on the cleavage plane, resulting in a ID wire (Figure 2a). We are able to apply a 
magnetic field either parallel or perpendicular to the wire (in the x or y directions indicated in Figure 2b). 

The basic properties of our wires have been or will be reported in other publications [2l| . Another body of work on 
quasi-one-dimensional wires in GaAs by Hamilton et al. [23L[2~3| focused on gate-defined quantum point contacts on a 
2DHG grown in the 311 direction. These studies found that, due to the SOI, the effective g-factor depends strongly 
on the particular subband studied as well as the direction of the wire and that of the magnetic field with respect 
to the crystal axes. They did not, however, observe signs of the SO gap reported below, perhaps due to the growth 
direction of their 2DHG or their nearly-symmetrical confinement potential. 
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FIG. 2: The device and measurement setup, a, Cross-section of our devices, which are fabricated by the cleaved-edge overgrowth 
method [2p| . l2ll ]. b, Measurement scheme for the one-dimensional hole wires. A section of the wire is isolated using a gate 
which depletes the two-dimensional hole gas (2DHG) just beneath it. Conductance is measured using ohmic contacts to the 
2DHG on either side of the wire and decreases in steps each of which corresponds to the depletion of a ID subband [21| 
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FIG. 3: The first two conductance steps of a quantum wire (blue traces) and their evolution in magnetic fields applied in two 
different directions (red traces), a, A field perpendicular to the wire splits the second step into two half-steps, the first marked 
with an arrow), b, A field parallel to the wire produces a half-step (arrow) and a dip (fat arrow), signifying the presence of 
a spin-orbit gap. On the first step (gate voltages larger than 3.5V) no dip is observed and half-steps, if they exist at all, are 
hardly discernible. 



In Figure 2b, applying a positive voltage to pre-fabricated gates on the top surface of the wafer, we deplete first the 
2DHG under the gate to isolate the ID wire and subsequently the subbands of the wire. The conductance, measured 
between ohmic contacts to the 2DHG on either side of the gate/wire, decreases in steps as the wire subbands are 
successively depleted (25, 26|. 

In Figure 3 — which shows our main experimental results — we focus on the first two conductance steps seen at 
zero magnetic field (blue traces) corresponding to the lowest two subbands in the wire, and their evolution in magnetic 
fields applied in two different directions, y and z (red traces). Concentrating first on the second, higher, conductance 
step, we see that a magnetic field in the y direction transforms this step into two half-height steps (Figure 3a), 
whereas a field in x produces both two half-steps and a dip (Figure 3b). In contrast, the first step appears unaffected 
by magnetic field in either direction. 

Let us try to understand these results in terms of the simple model presented in Figure 1. Assuming that the main 
contribution to the spin-orbit effect is due to structural inversion asymmetry, we expect B so to be perpendicular to 
the wire, in the y-z plane, as the structure is translationally invariant along the length of the wire. Thus, a magnetic 
field in the x direction (perpendicular to B so) should produce dips in the zero-field conductance steps (as in If), 



4 




FIG. 4: Predictions from our model of band structures and conductances traces, a, The dispersion relations of the two lowest 
quantum wire subbands with no applied magnetic field, b, The same bands in a magnetic field of 9T perpendicular to the wire, 
with g — 2. c, The same bands in a magnetic field of 9T parallel to the wire, with g = 2/9. d, Conductance traces calculated 
from a (blue) and b (red) with the experiment temperature of 300mK. Compare to Figure 3a. e, Conductance traces calculated 
from a (blue) and c (red) with the experiment temperature of 300mK. Compare to Figure 3b. 



whereas one in the y direction should split the each zero-field conductance step into two half-steps (as in le) and 
possibly also produce dips simultaneously. 

Our data bear more than a passing resemblance to the expectations from the simple model, yet the two are 
significantly different in their details. A more realistic model, described below, allows us to more fully understand our 
results. 

We use a four band Luttinger model [l3L and include confinement in two directions, taking into account the 
orientation of our wire in the GaAs crystal. This model considers the lowest four spinful subbands, with the higher 
subbands being ignored. The parameters in the model are: C, 71, 72, 73, Xi d y , d z , r y and r z . C and the ji are 
bulk (Al)GaAs band structure parameters which are well-established in the literature, x describes the leakage of the 
ID hole wavefunctions from the GaAs where they are nominally confined into the surrounding AlGaAs; however the 
model is not very sensitive to this parameter. d y and d z describe the confinement in the y and z directions, and r y 
and r z the strength of the SIA electric field in the y and z directions. We start by using an isotropic value of 2 for 
the g-factor and return to this point later. The full Hamiltonian and details of its derivation as well as values chosen 
for each parameter can be found in Supplementary Information. 

Figure 4a shows the lowest energy bands calculated with this Hamiltonian in the absence of magnetic field. Focusing 
first on the upper pair of bands, we see that applying a field in the y direction shifts them in energy (Figure 4b), 
producing a double step in the conductance trace (Figure 4d), which is seen in our data (Figure 3a). In contrast, a 
field in the x direction produces a gap in the band structure (Figure 4c) and a dip in the conductance trace (Figure 
4e), also seen in out data (Figure 3b). 

As for the lower pair of bands, the spin-splitting at B = is much smaller than in the upper bands. Additionally, 
as in our data, magnetic fields in both directions do not have much of an effect (Figures 3 and 4) Q . 

Clearly the coupling between the two lower bands due to the SOI, or a magnetic field, is suppressed. This is a 
consequence of the fact that these bands have been chosen in the model to be primarily of heavy-hole character — 
their carriers have spin ±3/2 — a choice motivated by our data and other physical considerations as explicated in 
the SOM. The leading order Rashba and Zeeman terms are linear order in spin-operators and thus can only couple 
states with AS — 1. As AS = 3 for the heavy holes, the two lower two subbands are not coupled by these terms. 
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FIG. 5: Reproducibility. Additional data, from a separate device, showing the dip on the second step due to the presence of 
the spin-orbit gap. 



Their coupling requires cubic order terms in the spin operators. Such terms, however, are small in all components 
of the Hamiltonian: Rashba and Zeeman couplings which are cubic in spin-operators do exist, but their coupling 



constants are suppressed [13|. The lowest order Dresselhaus term is cubic in spin-operators, but in GaAs it is again 



very small [13J and thus has almost no effect on the bands. 

Indeed, the Dresselhaus term has a minimal effect not only on the lower bands but on the upper pair of bands as 
well. Therefore, the SO effects seen in our samples come primarily from the asymmetry of the confining potential of 
the ID wire, i.e. the Rashba term. We find further that the electric field due to this potential is stronger in the z 
direction than in the y. 

We return now to the question of the g- factor. We find that, in order to produce features of the experimentally 
observed sizes at the applied 9T field, our model requires an effective g-factor of 2 in the y direction and an effective 
g-factor of 2/9 in the x direction. (Figure 4d and e.) Such an anisotropy in the effective g-factor is expected in the 
presence of SO coupling and is dependent on the strength and direction of confinement HI, 2^, 3, 3| . 

We are thus able to describe most of the important features of our data with the proposed model and to identify 
their physical causes. Further theoretical work is needed to understand one feature of our data which our model fails 
to capture: the occurrence of a double step in addition to a dip when the magnetic field is applied in the x direction. 

We also acquired data similar to those in Figure 3 at finely-spaced intermediate fields (Supplementary Information). 
These data show that the observed conductance dip and half-steps develop gradually with magnetic field, allowing us 
to rule out resonances associated with particular values of the magnetic field as the cause of the observed features. 
In addition, these features — and in particular the conductance dip — were observed in a separate device (Figure 5) 
and persist after thermal cycling; thus they are also not linked to particular configurations of the disorder potential 
in the device. 

In conclusion, we have observed for the first time the spin-orbit gap predicted for one-dimensional systems and we 
have developed a model which describes most of the detailed features of our data by taking into account the materials 
properties of (Al)GaAs. 

Several directions seem promising for the further exploration of these gaps. For this particular device geometry, 
even cleaner wires would allow the study of SO effects at finite bias voltage and in higher subbands. A tunnelling 
setup such as in Ref [1(3] could probe the dispersion relations directly; and a side gate as in Ref [^H would allow 
tuning the strength of the electric field which gives rise to the spin-orbit effect. 

For the purposes of spintronics applications, verification of spin transport, particularly at zero bias, would be of 
great interest. This could be achieved either through direct detection of the spin current or though the detection of 
spin accumulation at the two ends of the device. It would also be of practical interest to explore the possibility of 
producing such SO gaps in the lowest subband and without a magnetic field. While in our experiments the spin-orbit 
gap is induced by a magnetic field, this may not be necessary: the two spin bands could be mixed at k — in 
other ways, perhaps through the use of materials with magnetic order or through controlled doping with magnetic 
impurities. 

Finally, in other studies Coulomb interactions have been found to have marked effects in one-dimensional systems 
including CEO electron wires [HI, and further work is necessary to understand their role in these newly-developed 
systems. 
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SUPPLEMENTARY INFORMATION 
Derivation of the Hamiltonian and Choice of Parameters 

We begin with the four-band Luttinger model, which models the uppermost valence bands of GaAs This 
model, containing the heavy-hole (HH) and light-hole (LH) valence bands, is valid for large-gap III-V semiconductors 
with band edges at the T-point. For the bulk Hamiltonian we have 

h 2 x + 

k 

/-(7i+72)fcy -(71-272)^ 2V37 3 fc_fe y/3(ya& - i2j 3 k x k y ) ^ 

2V3 73 fc+fc 2 -(71 -72)Af -(71 + 272)^ Va(r/2& ~ i^akxky) 

V3(^ 2 K + i2^ 3 k x k y ) -(7i-7 2 )fc2_( 7l+ 2 72 )fc2 -2V3 73 fc_fc z 

V V3(j 2 K + i2 73 k x k y ) -2V3 73 fc+fc z -( 7l + 72 )fefj - (71 - 2 l2 )k 2 J 

Ck 

(1) 

k\ = k 2 x + k 2 (2) 



k = k 2 x - e y (3) 



k x ± iky (4) 



c fe = (c3/2k Ci/ 2 fc C_i/ 2 fc C_ 3 /2fc) (5) 

Here c a k destroys a fermion in the J z = a bulk band. As mentioned in the main article, we have found that the 
effects of bulk inversion asymmetry are very small and therefore inconsequential; we have therefore set the parameter 
C in Reference [l3| to zero to simplify the representation. The quantum well growth directions are the crystal [001] 
and [110] directions so we will rotate our axes around the [011] direction by an angle of -7r/4. From here on the x, y 
directions will be aligned along the [110] and [110] crystal directions respectively. 

In our device, holes are confined in the y and z directions, to create a ID wire along x. Therefore, ky and k\ may 
be replaced by their expectation values in the lowest subband: (k^) ~ (n/da) 2 where d a is the confinement width in 
the a direction. Ignoring structural inversion asymmetry (SIA) terms, the effective ID Hamiltonian is: 



h 2 ^ 



/ -(7i+72)^ -i{Vil2k 2 x -f3y) _0 \ 

-( 7l - 72 )fc2_A -i(V3 l2 k 2 x - (3 y ) 

i{VZl2k 2 x -P v ) _ _( 7l _ 72 ) fc 2_ A 

V *(V3 72 fc2 - (3 y ) _( 7l+72 )fc2 J 



ck (6) 
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with 



472 7T 



(7) 



Py = ^72^2" 



and the basis 



(8) 



( c 0+3/2fc c 0+l/2fc c 0-l/2fe c 0-3/2fe) 



(9) 



Here co CT fc destroys a fermion in the lowest confinement subband of the J = a bulk band at momentum k. We have 
chosen a pair of subbands arising from heavy-hole like states Co±3/2fc, and a pair of subbands arising from light-hole like 
states co±i/2fc- Although it is not clear a priori which bulk bands would give rise the four lowest subbands, this choice 
is motivated by the fact that, given Equation (7), the spacing between the chosen subbands could be comparable 
to the inter-subband spacing of the heavy-hole-like series of subbands. Thus, the two reasonable options seemed 
to be either that the two lowest sub-bands should be heavy-hole-like or that one of them should be light-hole-like. 
Our choice of the latter is subsequently validated by our data (see main article). The parameter (3 y arises from the 
confinement in the y direction [131 ]. 

We now return to the SIA (Rashba) terms. These may be deduced from symmetry, or more formally the method 
of invariants [l3j |. For spin-3/2 systems there are two allowed Rashba terms which are linear in momentum: 



Hsu 



ax ( 



k x E) ■ J 



U2 



J' 



(10) 



with J = (J x , J V ,J Z ) the spin-3/2 matrices and J = (Jj, Jy, Jf ). Empirically a-i <C a.\ so we shall ignore the second 
term For the confinement electric fields we have E z 7^ and E y 7^ and thus the ID Rashba term is 



ttR 



a x k x {E y J z - E z J y ) 
( 



^-r z k x 



V 



2 T*y kx 

^-r z k x 






^+r z k x 

— lfyk x ) 



(11) 



where r y = a\E, 



a\E z and <p± = 1 ± i. 



Finally, we consider terms added to the Hamiltonian to describe the effects of an applied external magnetic field. 
We shall ignore here orbital magnetic field effects because our magneto-conductance data does not reveal substantial 
shifts of the band edges with field, which are to be expected had orbital effects been important. As for the Zeeman 
terms, two are expected for the spin-3/2 systems: 

H z = g lf i B B ■ ~J + ,92MB B ■ f. (12) 
Here again gi -C gi and therefore the second term can be ignored [HI]. In matrix form, the remaining (first) term is: 



H z = g^B 



^4>+B + \B Z 
^+B + 



\ 



_ 

_ 



\ 



2 ^z 



^+B, 



(13) 



\B Z ) 



with g the Lande g-factor, fig the Bohr magneton, B± = B x ± iB y and cf>± = 1 ± i. 
Our full model Hamiltonian is: 



H — Hid + H ld + H z 



(14) 
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We have calculated the resultant dispersion relations for these four low-lying bands using a sensible set of parameters 
listed below and deduced the magneto-conductance traces expected from this model. Despite the various assumptions 
and approximations made, this model captures faithfully most of the features in our data semi-quantitatively (see 
main article). One feature in our data is not captured by this model; however, pointing to a particular assumption or 
approximation made here as the primary cause of this shortcoming is beyond the scope of this manuscript. 

We now turn to the values used for the various model parameters. We have used known values where possible and 
the reasoning for the choices made for the remaining parameters is detailed below. The values used are: 



TABLE I: Parameter values used in the model. 



Parameter 


Value used in model 


7i 


6.85 (5.68) 


72 


2.1 (1.63) 


73 


2.9 (2.35) 


X 


0.05 


dy 


60nm 


d z 


25nm 


r y 





r z 


23±2meV-nm 



The 7„ are crystal parameters and are well-established in the literature [131]. The numbers in brackets are for 
Alo.324Gao.676As calculated using a linear interpolation in the virtual crystal approximation, while those not in 
brackets are for GaAs. The actual numbers used in the model are 7„ = X1n,AlGaAs + (1 — x)ln,GaAs where x 
represents the leakage of the confined subband wavefunction into the AlGaAs (effectively modelling deviation from 
an infinite potential well). We estimated x from simulations of the wavefunction of the lowest bound state in the 
potential wells in both confinement directions using the program developed by Gregory Snider for this purpose (33[, 
exploring a rather large range of possible effective masses (0.2-0.65mo, where mo is the bare electron mass), guided 
by theoretical values give n in Reference [l3[ and the little that is experimentally known about the effective masses in 
these systems 34, 3^, 



361 ] . We obtained < x < 0-09 an d chose a value in the middle for the model. We note that x 
is in any case a parameter to which our results are rather insensitive. 

As noted above, C in Reference 13j, which has a value of 0.009271, has been set to zero as it is in any case too 
small to perceptibly affect the results of the model. 



The remaining four parameters dy , d z , r y 



are used as fitting parameters. They are grouped into two two- 



parameter sets which are treated differently; d y and d z , the confinement related parameters are constrained to rea- 
sonable values — not too different from an estimated extent of the wave function along these two directions, as 
explained below. We are uncertain about the values of d v and d z because while the thicknesses of the layers giving 
rise to the confinement are known to great accuracy, the actual extent of the wavefunctions and the precise relation of 
those to d y and d z are difficult to calculate. The other two parameters, r y and r z , the Rashba coefficients are simply 
tuned to match the model calculations to the data. The values found are in fact the first available estimates of these 
parameters. This is because most work on SO in the valence band of GaAs has focused on heavy holes, whereas these 
are the parameters relevant for the light- hole- like subbands. As noted in the main article, we find the Rashba field to 
be much stronger in the z direction than in the y direction. 

We have estimated the extent of the wavefunction in each direction rather crudely — by modelling two fictitious 
two dimensional hole gases separately, neglecting the joint effect of both confinements in the real device. The two 
2DHGs were simulated using nominal growth parameters in each direction separately. We found the extent of the 
wavefunction in z to be about 15nm, while in y, we obtain about 30-35nm. The values for d y and d z used in the 
model are within a factor of two of these figures. 

Finally, we note that our nonlinear transport measurements in the same device (to be published elsewhere) point to 
a subband spacing of about 3meV, which is slightly smaller than what the model predicts with the chosen parameters 
(~4meV, Figure 4a of the main article); we suspect that this mismatch is due to the mixing with nearby subbands 
(next few conductance steps), which are not considered in this model. 
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Figures a and b show data from Figure 3 in the main text, with additional traces taken at intermediate fields (every 
0.2T) showing that the features discussed in the text (arrows) develop gradually with field and are not spurious effects. 
Figure c shows the data in Figure b around the spin-orbit gap. 



